Is the Dark Disc contribution to Dark Matter Signals important ? 



Fu-Sin Ling 1 

1 Service de Physique Theorique, Universite Libre de Bruxelles, 
Boulevard du Triomphe - CP225, 1050 Bruxelles - Belgiqu^ 

Recent N-body simulations indicate that a thick disc of dark matter, co-rotating with the stellar 
disc, forms in a galactic halo after a merger at a redshift z < 2. The existence of such a dark disc 
component in the Milky Way could affect dramatically dark matter signals in direct and indirect 
detection. In this letter, we discuss the possible signal enhancement in connection with the char- 
acteristics of the local velocity distributions. We argue that the enhancement is rather mild, but 
some subtle effects may arise. In particular, the annual modulation observed by DAMA becomes 
less constrained by other direct detection experiments. 
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I. INTRODUCTION 

The problem of determining the local phase-space 
structure of Dark Matter (DM) in the solar neighbor- 
hood is rather difficult. The distance of the Sun to the 
galactic center is Rq = 8.0 kpc, which represents only 
a small fraction (around 3-4%) of the virial radius of 
the Milky Way galaxy. One would therefore expect most 
DM particles to have virialized, with only a very small 
fraction in cold flows, i.e. flows with negligible velocity 
dispersion [31 13] • This is indeed confirmed by high resolu- 
tion N-body simulations, which find only faint evidence 
for caustics [5]. 

However, the distribution of DM in galactic halos 
is not smooth. Hierarchical clustering in the ACDM 
paradigm leads to nested structures with numerous sub- 
halos Although recent DM only simulations [3 E] 
have reached a tremendous mass resolution, they are still 
far from resolving the smallest substructures that are ex- 
pected on theoretical grounds, with a mass in the range 
10~ n • • • 10~ 3 solar masses, depending on the decoupling 
temperature of the WIMP [9tiT2], 

When baryons enter the game, the presence of an al- 
ready well formed galactic disc at high redshift perturbs 
the merger process of the host galaxy with smaller satel- 
lites, as it causes a preferential drag in its direction [13] . 
As a result, the material from the accreted structure, 
including DM, stars and gas, settles around the galactic 
disc of the host halo, after being disrupted by tidal forces. 
In particular, the accreted DM forms a thick disc, often 
dubbed the dark disc, which co-rotates with the stellar 
disc. 

A consequence of this scenario is the formation of a 
thick disc of stars in the host galaxy, as the result of 
the heating of the pre-existing thin disc of stars during 
the merger, supplemented by the accreted stars. Many 
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morphological and dynamical properties of the accreted 
stars and DM should therefore evidence some correla- 
tion [TU US] ■ Using three disc galaxies from cosmological 
hydrodynamics simulations, that are intended to span a 
large range of merger histories, from very quiescent to 
very violent, Read et al. [I3J conclude that the existence 
of a thick stellar disk and a co-rotating dark disk in the 
ACDM paradigm is quite robust. 

The detailed quantitative properties of the final galaxy 
depend however on its merger history. In particular, the 
authors find that the dark disc can contribute ~ 0.25—1.5 
times the non-rotating halo density at the solar position. 
Also, the rotation lag of the rotating component is in 
the range — 150 km/s. A large value of the rotation 
lag, corresponding to a small halo circular speed, is only 
found for galaxies which had no significant merger after 
a redshift 2 = 2. 

The dramatic influence of the dark disk on DM signals 
has been shown in Ref. |16l for direct detection and in 
Ref. [E] for indirect detection. For direct detection, the 
low velocity of the dark disc with respect to the Earth 
enhances detection rates at low recoil energy. This could 
lead to an increase by up to a factor of 3 in the 5 — 20 keV 
recoil energy range, for DM masses bigger than 50 GeV. 
The annual modulation is also affected, with a reduction 
of the phase corresponding to maximal rate by up to 
three weeks. For neutrino indirect detection, the boost 
in signal is even more impressive. The enhanced capture 
rate in the Sun and in the Earth leads to a muon flux from 
the Sun (from the Earth) that can be one (three) orders 
of magnitude higher than in the Standard Maxwellian 
halo case. 

To arrive to their conclusions, the authors in Ref. [TS] 
and Ref. [17] assumed a lag velocity vi ag = 50 km/s, and 
supposed that the velocity distribution of the dark disc 
can be correctly described by an isotropic Maxwellian dis- 
tribution with a velocity dispersion a = 50km/s, neglect- 
ing anisotropics and deviations from Maxwell-Boltzmann 
statistics, and also taken to have the same value . They 
argue that this value of the velocity dispersion is sup- 
ported by the observed kinematics of the Milky Way thick 
disc stars, (an, a^, a z ) = (63,39,39) ±4 km/s [18], and 
is compatible with the results of numerical simulations. 
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The purpose of this letter is to discuss whether these 
assumptions are realistic or not, and what would the im- 
plications be for direct and indirect detection. 

II. DYNAMICAL PROPERTIES OF THE DARK 
DISC 

The dark disc affects the predictions of direct and in- 
direct detection signals essentially because it increases 
the number of dark matter particles that move slowly 
with respect to the Earth or the Sun. The most relevant 
parameters that determine the phenomenological impor- 
tance of the dark disc are therefore its local density po, 
its lag velocity vi ag with respect to the stellar disc, and 
its velocity dispersion a. 

One can estimate the velocity dispersion by applying 
the virial theorem. For a galaxy of the size of the Milky- 
Way, 

(v 2 ) = M ~ (200 km/s) 2 , (1) 

where W is the gravitational potential energy and M is 
the mass of the galaxy [19] . Therefore, if we assume a 
universal density profile for DM halos, a satellite halo 
with a mass M s around ten times smaller than the host 
halo mass Mh should have a velocity dispersion given by 

(v*) s ~ ^ (200 km/s) 2 , (2) 

which gives an average one dimensional velocity disper- 
sion around 50 km/s. This value is the mean velocity 
dispersion for the entire satellite halo, the question is 
to know how it is related to the velocity dispersion of 
the accreted DM at the Earth's location after merging of 
the satellite halo. The merging and disruption processes 
are likely to cause some heating. Also, the final veloc- 
ity dispersion will depend on whether the accreted DM 
originates rather from the inner or the outer parts of the 
satellite halo, as the average velocity of the particles drop 
with the distance from the center of the halo. 

In Ref. [13] , three cosmological hydrodynamics simula- 
tions of bright, disc dominated galaxies are used to pre- 
dict the phase-space distribution of dark matter in the 
solar neighborhood. The three galaxies have been cho- 
sen to span very different merger histories, from a very 
quiescent one after a redshift z = 2 to a violent one with 
a massive merger after redshift z = 1. This sample is 
meant to cover the possible extreme cases of the merger 
history of the Milky Way, and to bracket the characteris- 
tics of its dark disk. Without surprise, the contribution 
of the dark disc to the local dark matter density near the 
Sun increases when there are more accretion events. The 
authors find that pu — 0.25 . . . l.bpu, where pjj is the 
local density of the host halo. The lag velocity, between 
and 150 km/s, is also smaller for more massive mergers. 

Interestingly, on Fig. 4 of Ref. [13], the authors com- 
pare the distributions of accreted stars and dark matter. 



It appears that the accreted DM velocity dispersion is 
larger than 100 km/s for all three galaxies considered in 
their paper. More precisely, the quadratic average of the 
velocity dispersion in the directions R, </>, z amounts to 
a = 119.8, a = 108.1 and a = 135.7 km/s for the galaxies 
MW1, H204 and H258 respectively. For accreted stars, 
the average velocity dispersion is much smaller, with the 
largest value in the radial direction, a = 85.3, a = 83.5 
and a = 90.7 km/s for MW1, H204 and H258 respec- 
tively. As we can see, the velocity dispersion of the ac- 
creted DM is at least 25% higher than that of the accreted 
stars. This result is in agreement with the conclusions of 
Ref. [15] . where the authors use focused, high resolution 
simulations of accretion events to compare the morpho- 
logical and kinematical properties of the dark disc with 
those of the Milky Way's thick disc of stars. 

In Ref. [T5J, single accretion events at different orbital 
infall inclinations are studied with focused numerical sim- 
ulations. It appears that the velocity dispersion of ac- 
creted stars is substantially larger than the velocity dis- 
persion of the stars in the resultant thick disc of the host 
halo, except for very low inclination mergers. However, 
these planar events produce stellar discs that are much 
hotter than what is observed for the Milky Way [15] . 
Indeed, the thick disc of the Milky Way has a velocity 
dispersion (o\r, tr^, a z ) = (63 ± 6, 39 ± 4, 39 ± 4) km/s, it 
also has a moderate rotational lag vi ag — 51±5 km/s [18] , 
For higher inclination accretion events (9 ~ 60°), the re- 
sultant thick disc is less hot and compatible with obser- 
vations. However, in these cases, the accreted stars have 
a velocity dispersion 2 — 3 times larger than that of the 
thick disc, which results from the heating of the primary 
disc. Therefore, taking the observed values of the veloc- 
ity dispersions of the Milky Way thick disc as a way to 
characterize the dark disc velocity dispersion is unlikely 
to be realistic. Nevertheless, this provides a very conser- 
vative lower bound 1.25(<tr, ct^, a z ) ~ (89,50,50) km/s. 

More realistic values of the dark disc velocity disper- 
sion <jd are provided by full cosmological hydrodynam- 
ics simulations. We already mentioned that <jd exceeds 
100 km/s for all three galaxies considered in Ref. [13] . 
Another high resolution simulation with baryons has 
been analyzed in Ref. [50]. There, the velocity distri- 
bution along <f) for DM particles in a ring 7 < R < 9 kpc, 
\z\ < 1 kpc showed that the fraction of particles co- 
rotating with the stellar disk with a small lag velocity 
cannot be large. A fit with a double Gaussian distribu- 
tion gave the following possible parameters for the dark 
disc component: pp = 0.25pn, vi ag = 70 km/s, and 
<70 = 85 km/s. At this point, it is important to stress 
that the double Gaussian fit gives an acceptable descrip- 
tion of the simulation velocity distribution only in the 
direction <f). Indeed, if one infers isotropic velocity dis- 
persions do for the dark disc and oh for the halo from 
this fit, such extrapolation would lead to velocity distri- 
butions in the radial and z directions that are much more 
peaked than the actual ones. They would be leptokurtic, 
instead of platykurtic, as seen in the simulation [20] . 




Figure 1. (Left column) Normalized velocity distributions for different dark matter halos, in the galactic frame (top) and with 
respect to the Sun (bottom). Black dotted : standard Maxwellian halo, Blue thin dashed : slowly rotating Tsallis halo, Red 
dashed : halo with a strong dark disc, Green solid : halo with a mild dark disc and po/pH = 1/3, Green dashed : halo with 
a mild dark disc and pd/ph = 1/1- (Right column) Enhancement of the phase-space distribution compared to the standard 
Maxwellian halo, for the same halos, as a function of the velocity. See text (Section III) for the parameters of the Tsallis, 
strong and mild dark disc halo models. 



The failure of this extraction and extrapolation is re- 
lated to the fact that equilibrated dark matter haloes 
show systematic deviations from Maxwcll-Boltzmann 
statistics. The presence of a long-range gravitational 
force suggests that this kind of system is better described 
by the nonextensive Tsallis statistics [2T] 122] . Indeed, 



for the DM particles in a spherical shell around 8 kpc 
in Ref. [2D], the velocity distribution is very well fit 
by a Tsallis distribution with parameters q ~ 0.8 and 
Wo = 265 km/s. We recall that the Tsallis distribution is 
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defined as 



N(v ,q) 



1 -(!-«); 



9/(1-9) 



(3) 



where N(vo,q) is a normalization factor. For ring parti- 
cles with 7 < R < 9 kpc, |z| < 1 kpc, we can get a rea- 
sonable fit of the velocity distributions as a sum of two 
isotropic Tsallis distributions if we take qa = qD = 0.7, 
vff = 300 km/s, v$ = 200 km/s, and v iag = 70 km/s. 
With this decomposition, the dark disc velocity disper- 
sion is <td = 117 km/s, much larger than that inferred 
from the double Gaussian fit. The same discrepancy also 
appears in Ref. [13], as a double Gaussian fit of the v<p 
distributions (see Fig. 2 of Ref. [13]) gives Of) = 50.3, 
a D = 76.0, and a D = 87.9 km/s for galaxies MW1, H204, 
and H258 respectively, values that are more than 30% 
smaller than the average velocity dispersions of the ac- 
creted dark matter (see above). In other words, the dou- 
ble Gaussian fit of the distribution underestimates the 
dark disc velocity dispersion. 

We conclude that the velocity dispersion of the Milky 
Way's dark disc is likely to be larger than 100 km/s. 
As we have shown, this is indeed the case for all the 
galactic halos comparable to the Milky Way that are ex- 
tracted from full cosmological simulations with baryons. 
These halos have undergone several accretion events dur- 
ing their merger history, which might explain why the 
resultant dark disc velocity distribution is hotter than in 
some focused and controlled single accretion event simu- 
lations. 



III. DARK MATTER SIGNALS FROM THE 
DARK DISC 

The presence of a dark disc in the galaxy can lead to 
an enhancement of dark matter signals in both direct 
and indirect detection. Much of this enhancement stems 
from an excess, compared to the standard Maxwellian 
halo predictions, of particles with a small relative ve- 
locity with respect to the Sun (and the Earth). The 
contribution of the dark disc to dark matter signals has 
been calculated in Ref. [16 for direct detection, and in 
Ref. [17] for indirect detection. Both papers assumed a 
Maxwellian profile (isotropic) for both the static halo and 
the rotating dark disc, with parameters (vg , Vq, viag) = 
(220, 50^, 50) km/s. As discussed in Sec. jllj the Milky 
Way's dark disc is unlikely to have such a small velocity 
dispersion. In the sequel, we will refer to this parameter- 
ization as the strong dark disc model, for which we also 
assume Pd/ph = 1 and pn = 0.3 GeV/cm 3 . So, in this 
model, the host halo is the standard Maxwellian halo. 

We also consider a mild dark disc situation, which 
models well the velocity distributions extracted from 
Ref. [20]. We take isotropic Tsallis distributions for 
both the static halo and the rotating dark disc, with 
parameters qn = qD = 0.7 and (vQ,VQ,vi ag ) = 



(300,200,70) km/s. Here, we let the dark disc to halo 
density ratio have two values, pd/ph — 1/3 as in the 
simulation, or Pd/ph — T m order to compare with the 
strong dark disc model. For the mild dark disk case, the 
velocity dispersion is much larger arj = 117 km/s. We 
recall that such a large value is supported by hydrody- 
namics cosmological simulations [131 115) and kinematical 
properties of the Milky Way's thick disc [TS] . 

Finally, in order to separate the effects due to the 
adoption of Tsallis distributions from those of the dark 
disc, we also consider a Tsallis halo with parameters 
vtf = 290 km/s, q H = 0.7, p H = 0.4 GeV/cm 3 , which 
is globally co-rotating with the stellar disc with a speed 
l'h = 35 km/s, so that this model halo has the same 
average co-rotation velocity as the mild dark disc model 
halo with pd/ph = 1/3. 

The velocity distributions of all these models are shown 
on Fig. [l] in both the galactic frame and the Sun-based 
frame, with w = v~Vq and vq = (10, 225, 7) km/s is the 
velocity of the Sun in galactic coordinates [53] . The ro- 
tating dark disc leads to an increase of particles moving 
slowly with respect to the Sun. As we can see, the en- 
hancement is quite sensitive to the dark disc characteris- 
tics. While for the strong dark disc, the enhancement fac- 
tor in the limit w — > is about 50, for the mild dark disc, 
it hardly reaches 4.5 when the two dark discs have the 
same local density. With the realistic ratio puj Ph = 1/3, 
the enhancement factor is about 2. For the mild dark disc 
models, there is also a significant bump at higher veloci- 
ties, it is due to the shape of the Tsallis distributions with 
q < 1 which are flatter and fatter than the Maxwellian 
distributions. It is worth noticing that the bump is due 
to the static halo part, as it is not modified when the den- 
sity of the dark disc is varied. As a result, a Tsallis halo 
with the same global rotation whose velocity distribution 
in the galactic frame follows closely that of the mild dark 
disc halo model fails to produce a big bump at high ve- 
locity in the Sun-based distributions (see the respective 
positions of the blue thin dashed and green solid curves 
in top and bottom right panels of Fig. [IJ. However, a 
static Tsallis halo can nevertheless produce such a high 
velocity bump as seen in the mild dark disc halo model, 
although the absence of co-rotation of the dark matter 
halo is less realistic in the presence of a baryonic disc. 



A. Direct Detection 

Direct detection experiments are aimed at measuring 
the nuclear recoil of a WIMP-like DM particle with a 
fixed target. The differential event rate of nuclear re- 
coils as a function of the recoil energy Er is conveniently 
factored as 124 



dU 
~dE~ n 



Pdm 



do 



M DM dE 



R 



v(E R ,t) 



(4) 



where pdm is the local density of DM, Me>m is the 
WIMP mass, da/dEn encodes all the particle and nuclear 
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Figure 2. Direct detection - Elastic scenario : DAMA allowed regions in the elastic scenario compared with exclusion limits of 
CDMS and XENON, in the case of a standard Maxwellian halo (top left), a slowly co-rotating Tsallis halo (top right), a strong 
Maxwellian dark disc halo (bottom left), and a mild Tsallis dark disc halo (bottom right). DAMA contours are given at the 90, 
99 and 99.9% CL. Stars indicate local best-fit points. All other exclusion curves are at the 99.9% CL. 



physics factors, and r](Eji,t) is the mean inverse veloc- 
ity of incoming DM particles that can deposit a recoil 
energy Er. The time dependence of the velocity distri- 
bution is induced by the motion of the Earth around the 
Sun, which leads to a seasonal modulation of the event 
rate [251 

Here we consider both the elastic and the inelastic scat- 



tering scenarios. The cross-section factor is written as 



da 



M. 



N(7 o (Zf p + (A-Z)f n y 



•I n 



F 2 (E I 



(5) 



where Mm ~ AM n is the mass of the target nucleus 
with atomic number A and Z protons, fi n is the re- 
duced WIMP/nucleon mass, er" is the zero momentum 
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Figure 3. Direct detection - Inelastic scenario ; DAMA allowed regions in the inelastic scenario compared with exclusion limits 
of CDMS and XENON, in the case of a standard Maxwellian halo (top left), a slowly co-rotating Tsallis halo (top right), a 
strong Maxwellian dark disc halo (bottom left), and a mild Tsallis dark disc halo (bottom right). DAMA contours are given at 
the 90, 99 and 99.9% CL. Stars indicate local best-fit points. All other exclusion curves are at the 99.9% CL. 



WIMP-neutron effective cross-section, f p and /„ are the 
effective coherent couplings of the WIMP to the proton 
and the neutron, and F 2 (Er) is the nuclear form fac- 
tor which characterizes the loss of coherence for non zero 
momentum transfer. We use for the form factor the sim- 
ple parameterization given by Helm [27J 12H]- For elas- 
tic scattering, we consider scalar interactions for which 
f p ~ f n . For inelastic scattering, we take the weak inter- 



action Z boson exchange as a typical realization [29H31] , 
with f p /f n ~ -0.08 and a Q n = & F tfj2*. 

The dependence upon the DM velocity distribution is 
encoded in the quantity 

v(E R ,t)= I d 3 v f( ? , (6) 

Jv min \V -«©(*) I 

where v and u® (t) are the WIMP velocity and the Earth 



7 



velocity in the galactic frame. A minimum impact veloc- 
ity v m i n is needed to provoke a recoil inside the detector, 

with /i the WIMP-nucleus reduced mass, and S is the 
mass splitting between the dark matter candidate and 
the next-to-lightest dark state in the inelastic scenario. 
Eq. ([7]) shows that, in the case of elastic scattering, the 
dark disc contributes to low recoil energy events. It dom- 
inates over the halo signal below some value of the recoil 
energy which increases with the WIMP mass for a fixed 
target. As this analysis has been carried out in details in 
Ref. [JJ)] for a strong dark disc, we will not repeat it here. 
In the case of inelastic scattering, the enhancement of the 
phase-space at low velocities produced by the dark disc 
will contribute to the direct detection signal only for very 
small values of 6, typically 8 < 30 keV. The most interest- 
ing range of values is however higher, 50 < S < 150 keV, 
as the inelastic scattering scenario was first devised as a 
new way to accommodate the annual modulation signal 
observed by DAMA [32], [33] with exclusion limits set by 
other experiments that are compatible with null results. 

Fig. [2] (elastic) and Fig. [3] (inelastic) show the DAMA 
allowed regions compared with limits from CDMS- 
Ge [34H36], CDMS-Si [35l [37], XENON10 [38l [39] and 
CRESST II gDl ST]. We refer the reader to these experi- 
mental papers as well as to Refs. [201 US] f° r a description 
of each experiment, and for values of the exposures, en- 
ergy thresholds, efficiencies, etc. that were used. We just 
mention that a more conservative limit for XENON10 has 
been chosen, following the analysis of Ref. [32] . 

1. Elastic scenario 



dark disc setup and in the case of a globally slowly rotat- 
ing Tsallis halo is the same, although the exact location of 
the channeling region is slightly different. We find that 
the improvement can be related to how fast the high- 
velocity tail of the distribution drops compared to the 
Maxwellian halo. So, the channelling region of DAMA 
can be more than marginally consistent with the other 
experiments if a realistic DM halo, with a mild dark disc 
component and platykurtic velocity distributions, as sug- 
gested by recent cosmological simulations with baryons, 
is considered. 



2. Inelastic scenario 

For the strong dark disc halo model as for the stan- 
dard Maxwellian halo, exclusion limits from CDMS and 
CRESST leave us with the same two solutions to the 
DAMA puzzle (see Fig. [3]). Indeed, for S > 30 keV, the 
minimum velocity i^mn for a recoil is above 200 km/s, 
so that the low velocity enhancement for the strong dark 
disc does not contribute to the signal. 

Things dramatically change when more realistic distri- 
butions are considered, either a slowly co-rotating Tsallis 
halo or the model with a mild Tsallis dark disc. While 
the low mass solution, Mdm < 100 GeV completely dis- 
appears, the high mass solution with Mr>u — few TeV 
is strongly enhanced. Again, the improvement in com- 
patibility between DAMA and other null experiments 
is mainly due to the sharper drop of the high velocity 
tail in Tsallis distributions. However, we notice that the 
mild dark disc setup gives a solution with a significantly 
smaller mass, compared to the globally co-rotating Tsal- 
lis halo. 



In the case of a strong dark disc, the so-called 
channeling region for DAMA [33] |33], with Mdm ~ 
10 ... 20 GeV, is the same as with the standard 
Maxwellian halo. Indeed, channeled events on Iodine 
(A t = 127) dominate in this region [321 133H35] . so that 
Vmin ~ 400 km/s in the recoil range 2 < Er < 6 keV 
where the modulation is observed. Therefore, the com- 
patibility of DAMA and other exclusion limits is not im- 
proved by the dark disc. The dark disc starts to dominate 
for DM masses above ~ 30 GeV. Compared to the stan- 
dard halo case, the so-called quenched region is displaced 
towards smaller cross-sections. It is however still totally 
excluded by other null experiments limits. 

In the case of a mild dark disc, the quenched region is 
only slightly affected. On the contrary, the flatter profiles 
of Tsallis distributions affects the channeling region sig- 
nificantly. This region now clearly protrudes outside the 
XENON10 exclusion limit towards smaller DM masses. 
While the best-fit point of the region is still excluded, 
some part of the 90% confidence level (CL) region is 
not. The level of improvement in compatibility between 
DAMA and the other experiments observed for the mild 



B. Indirect Detection 

Observation of very high energy neutrinos from the 
Sun or from the center of the Earth would be a convinc- 
ing proof of the existence of annihilating DM. The flux 
of these neutrinos is dictated by the capture rate C of 
the celestial body, rather than by the annihilation cross- 
section. Indeed, the number of DM particles inside the 
body obeys the equation [53] 



dN 
~dt 



= C- C A N 2 



(8) 



where Ta = CaN 2 /2 is the annihilation rate. Its solution 
is given by 



T A = ^-tanh 2 (i/T) 



(9) 



for constant C and C'a, with the time constant t = 
(CCa)" 1/2 . Typically, r - 10 s yr for the Sun and 
t ~ 10 11 yr for the Earth. The age of the solar sys- 
tem is around 4.5 billion years, therefore t/r ^> 1 and 
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Figure 4. Indirect detection : Enhancement factor of the DM capture rate in the Sun (left) and in the Earth (right) compared 
to the prediction for a standard Maxwellian halo, as a function of the DM mass. Red dashed / halo with a strong dark disc, 
Green solid : halo with a mild dark disc and Pd/ph = 1/3, Green dashed : halo with a mild dark disc and Pd/ph = 1/1. 



T A = C/2 for the Sun while tjr < 1 and T A oc C 2 for 
the Earth. 

The capture rate per unit shell volume by an element 
with mass Mjv, number density n and zero-momentum 
cross-section Co with the DM, is given by [2H US] 

dC p DM f Vmax fw(«oo) _ , v , 1fV . 
dV Mum Jo v <x> 

Mat fQmax 

vQ- e Jv)=a n-^ F 2 (Q)dQ ,(11) 

Z f X JQmin 

where fw is the normalized DM velocity distribution in 
the Sun-frame, at 'infinity' (i.e. neglecting the distor- 
tions induced by the gravitational potential of the celes- 
tial body) , v and DM velocity and the value of 
escape velocity for the given shell, they are related to the 
velocity at infinity by v 2 = f 2 sc + u L: v max = \//3^v eS c is 
the maximum velocity for which a capture is possible, fi 
is the WIMP-nucleus reduced mass, Q m in = Mdmv^/2 
and Qmax = P+Mjjmv 2 /2 are minimal and maximal 
values of the recoil energy for capture to be effective, 
F(Q) is a nuclear form factor such that F(0) — 1, and 
p± = 4M DM M N /(M DM ± M N ) 2 . 

To calculate the capture rate in the Sun, we use the 
density profile extracted from the reference standard so- 
lar model of Ref . [SD] , a chemical composition as found in 
Ref. [24], and an exponential form factor which largely 
reduces capture on heavy elements like iron. The es- 
cape velocity for the Sun ranges from v™" x — 1395 km/s 
in the center to v™ 1 ™ — 618 km/s at the surface. For 



heavy WIMPs, the impact of the escape velocity varia- 
tion on the capture rate is very significant. To calculate 
the capture rate in the Earth, we take a two layers (core + 
mantle) approximation model for the internal structure, 
with a different chemical composition for each layer [51] . 
The escape velocity for the Earth is considerably smaller 
(v™£ x = 14.8 km/s and u™ n = 11.2 km/s), so that form 
factors can be safely neglected. Here we consider only 
scalar interactions. Our purpose is to compare the pre- 
dictions between different halos, rather than get a pre- 
cise evaluation of the signal for a given halo. Therefore, 
we calculate the capture enhancement for models with 
strong and mild dark discs compared to the capture rate 
for the standard Maxwellian halo (Fig. [3| . In this way, 
we avoid many difficulties that plague the precise calcula- 
tion of the neutrino flux, in particular that coming from 
the Earth [17] . 

For the Sun, it appears that the capture enhance- 
ment is moderate, even with a strong dark disc. Our 
result seems to be in conflict with the the Fig. 1 of 
Ref. [T7], where enhancements factors larger than 10 oc- 
cur for masses smaller than 100 GeV. The authors agree 
on the fact that the enhancement that can be read off 
their Fig. 1, which is for the muon flux at the Earth's 
surface, should originate from a similar boost value in 
the capture rate, as the annihilation rate of dark matter 
has reached a steady equilibrium inside the Sun. Their 
calculation is entirely numerical, and details about form 
factors, or cross-sections are not given, which makes the 
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comparison difficult. 

In this work, we have focused on the capture boost 
for scalar interactions. For a dark matter mass around 
10 GeV, it appears that the capture boost CB is well ap- 
proximated by the following zeroth order analytical limit 



CB 



w ) 



PSH I(fw,SB.) 



(12) 



where p and f\y are the local density and normalized 
velocity distribution for the dark matter halo considered, 
whereas the same quantities with a subscript SH refer to 
the Standard Maxwellian halo, and 



I(fw) = 



fw(«oo) 



(13) 



Indeed, for Mom — 10 GeV, it can be checked that 
the capture rate is dominated by collisions on Helium, 
Oxygen and Iron, these are characterized by a maximum 
velocity v max larger than 700 km/s, and a form factor 
reduction smaller than 25%. The values of the capture 



boost as calculated with Eq. ( 12 ) are 4.47 for the strong 



dark disc halo model, and 2.49 and 1.44 for the mild dark 
disc halo models with Pd/ph — 1 and 1/3 respectively, 
in very good agreement with the exact values in Fig. [4] 
The capture boost becomes significantly larger than this 
approximation only for dark matter masses much smaller 
than the hydrogen atomic mass, or much larger than the 
iron atomic mass, as the maximum velocity for a capture 
Vmax decreases to values where the low velocity strong 
enhancement due to the dark disc comes into play. 

For the Earth, the capture enhancement is larger for 
the strong dark disc, up to a factor 50, and shows dips lo- 
cated around the resonances with nuclei that are present 
with a large abundance, mainly Oxygen, Magnesium, Sil- 
icon, Iron and Nickel. Finally, we note that for both the 
Sun and the Earth, the enhancement factor is limited be- 
low around two for the mild dark disc halo. Therefore, if 



future N-body simulations confirm that realistic DM ha- 
los for the Milky Way have similar characteristics than 
our mild dark disc model, then no spectacular boost is to 
be expected in neutrino indirect searches from the Sun 
and the Earth. 



IV. CONCLUSIONS & PERSPECTIVES 

In this letter, we have investigated the possible charac- 
teristics of the dark disc, in connection with results from 
N-body simulations with baryons. We argue that its con- 
tribution to the local density is mild, and that its velocity 
dispersion is quite large. Consequences for direct and in- 
direct detection are studied. Signal enhancements are 
limited, but some subtle effects lead to an improvement 
of compatibility between DAMA and other null experi- 
ments. 

Both direct detection experiments and neutrino tele- 
scopes probe the local DM environment, they are also 
both sensitive to the DM elastic scattering cross-section. 
However, a modification of the DM phase-space distribu- 
tion affects differently these signals. If the DAMA modu- 
lation signal is taken seriously, both the DM mass and its 
spin-independent cross-section would be more or less de- 
termined. This analysis shows how neutrino telescopes 
provide complementary information, the indirect signal 
would further give indications on the local halo struc- 
ture. 
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